Bragg spectroscopy of a cigar shaped Bose condensate in optical lattices 
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We study properties of excited states of an array of weakly coupled quasi-two-dimensional Bose 
condensates by using the hydrodynamic theory. We calculate multibranch Bogoliubov-Bloch spec- 
trum and its corresponding eigenfunctions. The spectrum of the axial excited states and its eigen- 
functions strongly depends on the coupling among various discrete radial modes within a given 
symmetry. This mode coupling is due to the presence of radial trapping potential. The multibranch 
nature of the Bogoliubov-Bloch spectrum and its dependence on the mode-coupling can be realized 
by analyzing dynamic structure factor and momentum transferred to the system in Bragg spec- 
troscopy experiments. We also study dynamic structure factor and momentum transferred to the 
condensate due to the Bragg spectroscopy experiment. 

PACS numbers: 03.75.Lm,03.75.Kk,32.80.Lg 



I. INTRODUCTION 

The experimental reahzation of optical lattices [l| is 
stimulating new perspectives in the study of strongly cor- 
related cold atoms Bose- Einstein condensates (BEG) 
in optical lattices are promising realistic systems to study 
the supcrfluid properties of Bose gases [1, 0, @ . When 
height of the lattice potential is low, the BEG have co- 
herency over all lattice sites. However, when height of 
the lattice potential increases so that nearest neighbor 
site hopping becomes difficult for atoms and then co- 
herence of the BEG is destroyed. The atomic conden- 
sate placed in a three-dimensional (3D) optical lattice 
can be best described by the Bose- Hubbard model [1]. 
The Bose- Hubbard model has been realized and quantum 
phase transition from weakly interacting Bose superfluid 
to a strongly interacting Mott insulator state was indeed 
observed experimentally 0, Q. Similar quantum phase 
transition was also observed when a 3D BEG was placed 
in a one-dimensional (ID) optical lattice 3. The insu- 
lator phase in the ID lattice is called a squeezed state. 
Kramer et al. llQ] have found the mass renormalization 
in presence of the optical potential. The renormaliza- 
tion decreases the value of the axial excitation frequen- 
cies which have been verified experimentally There 
are several theoretical calculations for dynamic structure 
factor of quasi-one-dimensional (ID) Bose gas placed in 
ID optical lattices [H, [H, Recently, Zobay and 

Rosenkranz [3] have studied effect of weak ID periodic 
potential on the thermodynamic properties (critical tem- 
perature, critical number and critical chemical potential) 
of a homogeneous interacting Bose gases by using mean- 
field and renormalization group method. Sound veloc- 
ity in an interacting Bose system in presence of peri- 
odic jDptical lattices has been studied by various authors 

0,1121 HsHi^. 

A stack of weakly coupled quasi-two-dimensional 
(quasi-2D) condensates can be created by applying a 
relatively strong ID optical lattices to an ordinary 3D 
cigar shaped condensate. The finite transverse size of 



the condensates produces a discreteness of radial spec- 
trum. The short wavelength axial phonons with differ- 
ent number of discrete radial modes of a cigar shaped 
condensates placed in optical lattices give rise to multi- 
branch BogoUubov-Bloch spectrum (MBBS) This 
is similar to multibranch Bogoliubov spectrum (MBS) of 
a cigar shaped BEG 

Bragg spectroscopy of a trapped BEG has proven to be 
an important tool for probing many bulk properties such 
as dynamic structure factor [23| , verification of multi- 
branch Bogoliubov excitation spectrum [24| , momentum 
distribution and correlation functions of a phase fluctu- 
ating quasi-lD Bose gases [1^ and even the velocity 
field of a vortex lattice within a condensate [131 . 

Recently, we have studied the MBBS by including cou- 
plings among all low energy modes in the same angular 
momentum sector by using hydrodynamic theory [20i |. 
This mode coupling is due to the presence of inhomoge- 
neous radial density. We should mention that the MBBS 
for phonon and breathing modes have been studied by in- 
cluding couplings between phonon and breathing modes 
only by using variational method [1^. The mode cou- 
plings among all low energy modes must be considered 
in order to calculate correct eigen spectrum and the cor- 
responding eigen functions [20|. The MBBS softens as 
we increase strength of the optical potential. It indicates 
a transition from the superfluid state to the Mott insu- 
lating state along the optical lattice but the superfluid 
state remains in each quasi-2D condensates. 

To understand effect of the mode coupling on the spec- 
trum and the softening of these modes as we increase the 
laser intensity, we need to measure those modes. There 
is no theoretical study on how to measure these modes 
in the superfluid regime. The MBBS can be measured 
by using the Bragg scattering experiments via dynamic 
structure factor measurements. In this work, we start 
from the beginning systematically by presenting results 
of the MBBS and the corresponding wave functions of 
the density fluctuations. Then we calculate the dynamic 
structure factor and momentum transferred to the sys- 
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tern which will give us information about the spectrum. 

This paper is organized as follows. In Sec. II, we con- 
sider a stack of weakly coupled quasi-two-dimensional 
Bose condensates. Using the discrctized hydrodynamic 
theory, we calculate the multibranch Bogoliubov-Bloch 
spectrum and its corresponding eigenfunctions by includ- 
ing the mode coupling within a given symmetry. In Sec. 
Ill, we study the dynamic structure factor. In Sec. IV, we 
analyze the momentum transferred to the condensates by 
Bragg pulses. We also calculate time duration needed for 
the Bragg potential to observe the multibranch nature of 
the spectrum. We give a brief summary and conclusions 
in Sec. IV. 



II. MBBS OF A STACK OF BOSE 
CONDENSATES 

We consider that bosonic atoms, at T = 0, are trapped 
by a cigar shaped harmonic potential and a stationary 
ID optical potential modulated along the z axis. The 
Gross-Pitaevskii energy functional can be written as 

= J dV^Hr,z)[-^\/' + VUr,z) 

+ ||^(r,z)|2 + Kp(zMr,z). (1) 

Here, Vhoir, z) = (M/2)(w^r^ + uj'^z'^) is the harmonic 
trap potential and Vop{z) — sErSm'^{qoz) is the opti- 
cal potential modulated along the z axis where Er — 
?i^g§/2M is the recoil energy, s is the dimensionless pa- 
rameter determining the laser intensity and qq is the wave 
vector of the laser beam. Also, g = Airah^ /M is the 
strength of the two-body interaction energy, where a is 
the two-body scattering length. We also assumed that 
ujr >> uJz so that it makes a long cigar shaped trap. The 
minima of the optical potential are located at the points 

= Ji'/'Zo = jd, where d = tt/qq is the lattice size along 
the z axis. Around these minima, the optical poten- 
tial can be approximated as Vop{z) ~ {M/2)lu'^{z ~ Zj)'^, 
where the layer trap frequency is ujg — ^fiq^/M. In the 
usual experiments, the well trap frequency is larger than 
the axial harmonic frequency, lOs >> ujz- For a typical 
experiment [Tl| , the system parameters are cj^ — Stt x 90 
Hz, ^2 ~ 27r X 9 Hz, and A = (Zo/Stt = 820 nm. When 
LUs » ujz and the number of quasi-2D layers is quite 
large, we can safely ignore the harmonic potential along 
the z axis. Moreover, we also consider the atom number 
in each quasi-2D condensates to be constant and equal 
at each layer. The number of atoms in each layer for a 
typical experiment is '--^ 10"*. 

The strong laser intensity will give rise to a stack of 
several quasi-two-dimensional condensates. Because of 
quantum tunneling, the overlap between the wave func- 
tions of two consecutive layers can be sufficient to en- 
sure full coherence. If the tunneling is too small, the 
strong phase fluctuations will destroy axial phase coher- 
ence among the layers and lead to a new strongly corre- 



lated quantum state, namely Mott insulator state along 
the z axis and superfluid state in each quasi-2D conden- 
sate. One can give a rough estimate of the threshold laser 
intensity for the exotic quantum phase transition to be 
quite large p9| . 

In the presence of coherence among the layers it is nat- 
ural to make an ansatz for the condensate wave function 



as 



ip{x,y,z) ^^'ipj{x,y)f{z - Zj). 



(2) 



Here, Tpj{x,y) is the wave function of the quasi-two- 
dimensional condensate at a site j and f{z — Zj) is a 
localized function at the j-th site. The localized function 
can be written as 

Substituting the above ansatz into the energy functional 
and considering only the nearest-neighbor interactions, 
one can get the following energy functional: 

^0 ^ E / dxdy[~^i:]S/li:j +VUrMj\^] 



52D 



^ / dxdy^]^]Tpjil)j 



J ^ I dxdy[il;]^^i}j+'tl;]i;j+s\. 



(4) 



J,<5=±1 



where J is the strength of the Josephson coupling be- 
tween adjacent layers which is given as 



J 



dzf{z)[-—Wl + V,p{z)]f{z + d) 
V 2A 



(5) 



where = yJfi/Mujr- Also, the strength of the effec- 
tive on-site interaction energy is g2D — 9 j dz\fQ{z)\'^ ~ 
AyJ 11 / 2{h'^ / M){a / as) , where as = \J1iJMlo's. Equation 
(III) shows that each layer j is coupled with the nearest- 
neighbor layers j ± 1 through the tunneling energy J. 
The axial dimension appears through the Josephson cou- 
pling between two adjacent layers. The Hamiltonian cor- 
responding to the above energy functional is similar to 
an effective ID Bose-Hubbard Hamiltonian in which each 
lattice site is replaced by a layer with inhomogeneous ra- 
dial density. 

Using phase-density representation of the bosonic field 
operator as V-'j = y^e'^^ and following Ref. [g^l, equa- 
tions of motion for the density and phase fluctuations can 
be written as 



2 7 

+ — no(r)[250, - 



(6) 



3 



and 



Using Eqs. (O, ©, ([8]), and dH), we get 



hS9j = —g2DSnj — 



2no{r) 



[2Snj — Srij^i — Srij-f-i]. (7) 



Here, " " represents the time derivative. In equilibrium, 
the condensate density at each layer within the Thomas- 
Fermi (TF) approximation is no(r) ~ [/ig — Vho('')]/52D, 
where we have neglected the effect of the tunneling energy 
J since it is very small in the deep optical lattice regime. 

Also, fiQ — fiLJr\J ^J&/n{Na/as) is the chemical poten- 
tial at each layer, where N is the number of atoms at 
each layer. Note that the second term (e(r) = J/2no{r)) 
of the right hand side of Eq. ([7]) is proportional to the 
small parameter J and inversely proportional to the large 
parameter no{r = 0) = /io/5'2D- In our previous publi- 
cation [l^l, we have neglected this small term. Within 
the TF regime, this term is really negligible. However, to 
present more accurate result we keep this term by approx- 
imating safely as J/2no{r) ^ J/2no{r = 0) = J(72d/mo- 
One may raise a objection that this approximation is not 
valid because the ratio e(r) is large as we go away from 
the trap center and it diverges at the boundary of the 
system. This divergence at the boundary is due to an 
artifact of the TF approximation to the density profile. 
The main idea of the linear response theory is to satisfy 
the following relation: 6n(r)/nQ{r) << 1. In our case, 
the term Jdn/2no{r) is always small within the linear 
response theory. 

The density and phase fluctuations can be written in 
a canonical form as 



and 



,(r,i) = ^[Sa.fe^„,fe(r)S„^fee'(^"'='*-""W*) + h.c.]. (9) 



Here, 5 is a destruction operator of a quasiparticle and a 
is a set of two quantum numbers: radial quantum num- 
ber, rir and the angular quantum number, m. Also, k 
is Bloch wave vector of the excitations. The Bloch wave 
vector p which is associated with the velocity of the con- 
densate in the optical lattice is set to zero. The density 
and phase fluctuations satisfy the following equal-time 
commutator relation: [6h(r)^ d9{r')] = i5{r — r'). One 
can easily show that 



' huJa{k) 



' 921) Bp 
2huOa{k) 



(10) 



where Bp — 1 + (2 J//io) sin^(fcc?/2) is a dimensionless pa- 
rameter. The parameter — 1 within the TF regime 
as well as in the tight binding regime. We are assum- 
ing that 'ipa,k{T) satisfies the orthonormal conditions: 

I d^r^a,k(^)^°"Ar) = ^aa' and Ea C.fe = 
Sir-r'). 



- ujl{k)tlja,k = So^Vr • [no(r)Vr'(/'a,fe] 

- Bo^^noir)sm\kd/2)^P^^k.in) 

For fc = 0, the solutions are known exactly and ana- 
lytically [30| . The energy spectrum and the normalized 
eigen functions, respectively, are given as = cj^(fc = 
0) = ujr[\i^\ + 2nr{nr + \m\ + 1)] and 

^a(r,0) = ^a,fe=o(r,0) ~ fl'"lp,(l'"l'0)(l - 2f2)e™'^. 

Here, P^'^\x) is the Jacobi polynomial of order n and 
(f) is the polar angle. The radius of each condensate layer 
is i?o = 2fj,o/Mu!^ and f = r/Rp is the dimensionless 
variable. 

The solution of Eq. pT|) can be obtained for arbitrary 
value of k by numerical diagonalization. For A: 7^ 0, we 
can expand the density fluctuations as 



(12) 



Substituting the above expansion into Eq. lTTI) . we ob- 
tain. 



= [Cjl{k) - Bp[\m\ + 2nr{nr + \m\ + \)]]ca 



BpBi sin^(/cd/2) ^ M^a'Ca' 



(13) 



where uJa — iOal^r and the dimensionless parameter B\ 
is defined as B\ = SJflp. Here, J = J/TiLJr and fip — 
fip/fiujr. The matrix element Maa' is given by 

AW= j (ffra,{r,4>){l-f'')i'c.{r,4>)- (14) 

The above eigenvalue problem is block diagonal with no 
overlap between the subspaces of different angular mo- 
mentum, so that the solutions to Eq. ([13]) can be obtained 
separately in each angular momentum subspace. We 
can obtain all low energy multibranch Bogoliubov-Bloch 
spectrum and the corresponding eigenfunctions from Eq. 
p3)) . Equations p3)) and (fT4)l show that the spectrum 
depends on average over the radial coordinate and the 
coupling among the modes within a given angular mo- 
mentum symmetry for any finite value of k. Particularly, 
the couplings among all other modes are important for 
large values of kd and Bi . It is interesting to note that 
the curvature of a spectrum strongly depends on the pa- 
rameter J5i and weakly depends on the parameter Bp. 
The tunneling energy (J) decreases as we increase the 
strength of the optical potential and therefore the modes 
become soft. 

When the radial trapping potential is absent, we obtain 
the usual Bogoliubov-Bloch spectrum from Eq. (|13p and 
it is given by 

^lomo{k) = 8 J/io sm^{kd/2) -f 16^ sin^(fcd/2). (15) 
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In the long- w avelengt h limit, the sound velocity is given 
by Chomo = \/ I^qIM* , where M* = Ti' j^jS is the effec- 
tive mass of the atoms in the optical potential. 

Now we consider inhomogeneous radial density. In the 
limit of long wavelength, the = 0, m = mode is 
phonon-like with a sound velocity cq = \J ^q/2M* . This 
sound velocity exactly matches with the result obtained 
in Ref. fTo*! and is similar to the result obtained without 
optical potential 21]. This sound velocity is smaller by 
a factor of \/2 with respect to the sound velocity for 
homogeneous condensate placed in ID optical lattices. 
This is due to the effect of the average over the radial 
variable which can be seen from Eqs. and (fH|) . 

In Fig. 1, we show few low-energy multibranch 
Bogoliubov-Bloch spectrum in the m = sector as a 
function of kd by solving the matrix Eq. (jl3p for given 
values of J = Q.lfwor and /ig = ?>()fwJr- Hereafter, we 
will also use these parameters for other figures. The tun- 
neling energy J ~ O.lhujr corresponds to the strength of 
the laser intensity s ^ 13. Within this parameter regime, 
atoms can tunnel from one layer to the adjacent layer and 
the whole system lies within the superfluid regime [29| . 
The lowest branch corresponds to the Bogoliubov-Bloch 
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FIG. 1: (Color online) Plots of the low-energy Bogoliubov- 
Bloch modes in the m — sector. Here, J — O.lhur and 

axial mode with no radial nodes. This mode has the 
usual form like to^k) = c^k at low momenta. The second 
branch corresponds to one radial node and starts at 2ujr 
for A: = 0. The breathing mode has the free-particle dis- 
persion relation and it can be written in terms of the ef- 
fective mass {ml) of this mode as W2(fc) = 2ujr+hk^ /2ml. 

The eigenf unctions corresponding to the low-energy 
modes are shown in Fig. 2. The eigenfunctions also sat- 
isfy the symmetry rule: ipn^.kir) = V'rv,27r-fe('')- 
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FIG. 2: (Color online) Plots of the eigenfunctions of the low- 
energy modes for J = O.lhur and fio = bOhuir- 

III. DYNAMIC STRUCTURE FACTOR 

The dynamic structure factor is the Fourier transfor- 
mation of density-density correlation functions. The dy- 
namic structure factor for the layers of quasi-2D conden- 
sates can be written as 

Siq,u;) ^ Yl j d^rd^r'dte"^'^'^^-^'^e'^* 

X < 5n]{r,t)5hj,{r' ,0) > . (16) 

Here, we have used the fact that q is perpendicular to the 
radial plane. At zero temperature, it can be re-written 
as 

S{q,i^) ^YSa{q)S{uj - UJa{q)), (17) 

a 

where weight factor Sa{q) is given by 

So^iq) = |A„,,n^„(g)p. (18) 

Here, ipa{q) = / <Pr^a,q{i')- The weight factors are 
shown in Fig. 3 as a function of kd. The weight factors 
Sa (q) determine how many modes are excited for a given 
value of q. For example, only rir — 1 and 2 modes are 
excited when qd = 1.0 and other modes are vanishingly 
small. The weight factors also satisfy the lattice symme- 
try, i.e., Sa{q) = Sa{2'K — q). In Fig. 4, we also plot 
the dynamic structure factor as a function of the Bragg 
frequency w for two values of the Bragg momenta q with 
fixed J and /iq. We find multiple peaks in the dynamic 
structure factor for a given value of q. The location of 
the peaks correspond to the excitation energy for a given 
value of q with fixed J and hq. These peaks show the 
multibranch nature of the Bogoliubov-Bloch modes. 

IV. BRAGG SPECTROSCOPY 

The behavior of these multiple peaks in the dynamic 
structure factor can be resolved in a two-photon Bragg 
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FIG. 3: (Color online) Plots of the weight factor for J 
O.lhujr and /^o ~ 507iaJr 
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FIG. 4; (Color online) Plots of the dynamic structure factor 
for kd = n/2 (solid) and kd = n (dashed) with fixed J = 
Q.lhuJr and no = SOfiWr- 



spectroscopy, as shown by Steinhauer et al. [2J]. The 
Bragg scattering experiments can be done by applying 
an additional moving optical potential in the form of 
Vsit) = Vbcos(gz — ujt). Here, Vb is the intensity of 
the Bragg pulses. The Bragg potential Vb (t) is indepen- 
dent from the static lattice potential Vop and it is also 
much weaker than the lattice potential. We also assume 
that the excitations are confined within the lowest band, 
therefore, uj must be less than the energy gap between 
the first and second bands. In the two-photon Bragg 
spectroscopy, the dynamic structure factor can not be 
measured directly. Actually, the observable in the Bragg 
scattering experiments is the momentum transferred to 
the condensate which is related to the dynamic struc- 
ture factor and reflects the behavior of the quasiparticle 
energy spectrum. The populations in the quasiparticle 
states can be controlled by using the two-photon Bragg 
pulses. When the superfluid is irradiated by an exter- 
nal moving Bragg potential Vsit) the excited states are 
populated by the quasiparticle with energy Hlu and the 
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FIG. 5: (Color online) Plots of the P^f) for J = Q.lUujr and 
/io = 50ftajr for various times t = 0.5rr(solid), t — l.OTr 
(dashed) and t = 2.0Tr (dot-dashed). 



momentum hq, depending on the value of q and to of 
the Bragg potential VB{t). Due to the axial symmetry, 
the modes having only zero angular momentum can be 
excited in the Bragg scattering experiments. Since the 
Bragg potential is much weaker than the lattice poten- 
tial, we can treat the scattering process with linear re- 
sponse theory. When the system is subjected to a time- 
dependent Bragg pulses, the additional interaction term 
appears in the total energy functional, which is given by, 

Hi{t)= dVi:''{r,z,t)[VoCOs{qz ~ ujt)]ip{r,z,t). (19) 



The Bragg potential can be approximated as VB{t) — 
Vb cos{jqd — ujt) since the atoms scatter off from the con- 
densate by the Bragg pulses from the j-th layer. The 
above interaction energy can be re-written as 

Hi{t) =J2[ tWo coa{jqd - ojt)i>j{r, t). 

Following the Rcfs. [Sll, [s^l, the momentum transfer to 
the Bose system placed in the optical lattices due the 
Bragg potential can be calculated analytically and it is 
given by 

P.(i) = < 6l,fe(t)6„,fe(t) 

a,k ^ 

X J2hqS^iq)iF^[u;_t]~Fo,[u;+t]), (20) 
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where ha.k{t) is the time-evolution of the quasiparticle 
operator of energy huja{k) and 



sm[{uja{q) ± Uj)t/2] 
K(g)±w)/2 



(21) 



In Fig. 5, we plot Pz{t) as a function of the Bragg 
frequency for two values of the Bragg momenta with var- 
ious time duration. For positive uj and a given q such that 
Saiq) is maximum, the momentum transferred Pz{t) is 
resonant at the frequencies io = uja{q)- The width of the 
each peak goes like 27r/<. Figure 5 shows that shape of the 
Pz {t) strongly depends on the time duration of the Bragg 
pulses. When t — 0.5Tr, the Pz{t) is a smooth curve for 
both values of the Bragg momenta. Here, = 27r/wr 
is the radial trapping period. When t — l.OTr, there 
is a clear evidence of few small peaks developed in the 
Pzit). When t = 2.0Tr, the multiple peaks in the Pz{t) 
appears sharply. The location of the peaks in Fig. 5 for 
t > l.OTr are exactly same as in Fig. 4. It implies that 
Pz{t) ~ S{q,iu) for a long duration of the Bragg pulses. 
Therefore, in order to resolve the different peaks, the du- 
ration of the Bragg pulses should be at least of the order 
of radial trapping period T^. By measuring the Pz at dif- 
ferent Bragg momenta, it is possible to get information 
about the MBBS. 



tides with various discrete radial nodes of an array of 
weakly coupled quasi-two dimensional Bose condensates. 
Our discretized hydrodynamic description enables us to 
produce correctly all low-energy MBBS and the corre- 
sponding eigenfunctions by including the mode couplings 
among different modes within the same angular momen- 
tum sector. We have also calculated the dynamic struc- 
ture factor and the momentum transferred to the system 
by the Bragg potential. In order to resolve the multiple 
peaks in the dynamic structure factor, the time duration 
of the Bragg pulses must be at least of the order of the ra- 
dial trapping period Tr- The MBBS can be measured by 
measuring Pz{t) for different values of Bragg momenta, 
similar to the experiment By measuring the MBBS, 
one may confirm the effect of the mode coupling on the 
spectrum. 
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